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A  METHOD  TO  COMPUTE  THE  FORCE 
SIGNATURE  OF  A  BODY  IMPACTING  ON  A  LINEAR 
ELASTIC  STRUCTURE  USING  FOURIER  ANALYSIS 

INTRODUCTION 

NRL  has  historically  been  involved  in  defining  shock  design  inputs  for  equipment  aboard  subma¬ 
rines.  A  body,  such  as  a  torpedo,  impacting  on  a  structure,  such  as  a  submarine  hull,  is  one  subset  of 
this  general  class.  This  report  presents  a  general  method  of  determining  the  force  signature  of  a  body 
impacting  on  any  linear  elastic  structure. 

Consider  equipment  which  has  been  attached  to  a  structure,  for  example  a  submarine  hull,  that  is 
struck  by  a  body  such  as  an  inert  torpedo  at  one  of  the  frames.  The  goal  is  to  determine  the  force 
exerted  by  the  body  impacting  on  the  structure  by  using  the  measured  response  at  various  gages  on  the 
equipment.  In  addition  the  impulse  response  at  various  gage  locations  is  needed.  The  response  to 
impulse  may  be  obtained  by  using  a  standard  computer  structural  analysis  code  such  as  NASTRAN. 

ANALYSIS 

Consider  responses  recorded  by  gages  at  several  points  on  the  equipment.  For  a  linear  elastic 
structure  the  response  at  point  P  due  to  a  force  at  point  K  is  the  convolution  of  the  force  at  point  K 
and  the  response  at  point  /*due  to  a  unit  impulse  at  point  K.  This  is  expressed  as: 

W')- Fk(T)  I„kU-  T)dT  (1) 

or: 

Rr-FK'  /,*.  (2) 

The  convolution  of  two  transforms  in  time  domain  is  the  inverse  transform  of  the  product  in  frequency 
domain.  Thus 

*,(«)  -  F*(«)  In t(»)  (3) 

where  FK(tt),  faitt)  are  the  Fourier  transforms  of  Rp(t),  FK(t),  and  /#*(/).  Rp(m)  for 

example  is  defined  by: 

/?,(«)  - RpiOe'1*'  dt  (4) 

and  its  inverse  transform  by: 

RpU)  -  *,(•»)*'*"  dm.  (5) 

In  order  to  make  use  of  a  very  accurate  numerical  method  to  compute  Fourier  sine  and  cosine 
transforms  derived  in  a  previous  report  U],  it  is  necessary  to  express  the  transforms  in  Eq.  (3)  in  terms 
of  sine  and  cosine  transforms  and  find  the  associated  inverse  transforms.  Consider  a  general  ftifiction 
f(p)  which  is  real  and  satisfies  the  Dirichlet  conditions  and  the  following  integral  exists 

/_7  l/Wlm  «> 
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Writing  the  Fourier  transform  of  fip)  as: 

Finn)  -  f(p)e~lmf  dp  -  J*_w  f(p)  cos  up  dp  -  i  f  ^  f(p)  sin  up  dp 
and  defining 

F(u)  *■  Fc(u)  —  i  Fsiu). 

Writing  the  inverse  Fourier  transform  of  F(u)  as: 


(7) 


(8) 


fit)  “  -jj1  J  ^  F(u)elm,du  *  (Fc  —  /^)  (cos  ut  +  i  sin  «/)  du  (9) 

then  rewriting 

/(r)  “  ^  J  Fc  cos  »r  </«  +  -j—  F$  sin  ut  du 

i  ii-fdi  f 

*  2t  ^  cos  «*f  sin  (10) 

Since  Fs(u)  —  -  Fsi- «*),  Fs  is  an  odd  function,  and  since  cos  ut  is  an  even  function  the  first  integrand 
on  the  second  line  of  Eq.  (10)  is  odd;  because  the  integral  is  over  the  symmetric  limits  to  -H»  it 
vanishes.  The  second  integrand  on  the  second  line  is  also  odd  and  therefore  both  imaginary  integrals 
vanish.  The  integrands  on  the  top  line  of  Eq.  00)  are  even  and  consequently: 

i  f  .  ..  .  ,  .  1 

ir 


fit) 


1  Cm  1  /*“ 

—  J  Fc  cos  u t  du  +  — J  Fssiautdu. 

*FF  *  0  w  (J 


(11) 


Since  for  a  physical  system  fip)  —  0  for  p  <  0  Fc(u)  and  Fs(u)  become 

Fc(o»)  «■  J0  fip)  cos  up  dp  (12) 

Fsiu)  -  fQ  fip)  mump  dp  (13) 

where  Fc  trnd  Fs  are  now  the  familiar  Fourier  cosine  and  sine  transforms.  If  fip)  is  an  even  function, 
from  Eqs.  (7)  and  (9): 


1  #•+*» 

fit)  -  —  J  m  fip)  cos  up  cos  c tit  dp  du 


2 

fcit)  •  ■—  Jo  Fciu)  cos  ut  du 


which  may  be  written  as: 

if 

it  J  o 

and  is  the  inverse  Fourier  cosine  transform.  This  transform  and  the  Fourier  cosine  transform: 
Fciu)  m  f0  fcip )  cos  up  & 

make  up  a  convenient  transform  pair.  If  fip)  is  an  odd  function  from  Eqs.  (7)  and  (9) 


fit)  «•  J_m  fip)  sin  up  sin  ut  dp  du 


which  may  be  written  as: 

2  rm 


(14) 


(15) 


(16) 


(17) 


(18) 


2  t 

f$(t)  —  —  J  Fgiu)  sin  ut  du 

and  is  the  inverse  Fourier  sine  transform.  This  transform  and  the  Fourier  sine  transform: 

Fsiu)  -  /0  fsip)  tin  up  dp  (1W 

make  up  a  transform  pair.  Returning  to  Eq.  (3)  and  writing  in  terms  of  sine  and  cosine  transforms 

(. He  -  /  *s)p  -  iFc-i  Fs)k  t/c  -  Ksl*.  (»> 
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2 


Dropping  the  subscripts  for  convenience  when  dealing  with  oaiy  one  response  point  in  this  deriva¬ 
tion,  yields,  when  setting  real  and  imaginary  parts  equal:  0 

Rc  m  FCIC  ~  Fg/g  (21) 

and 

Rs  m  FqIs  +  fj/c!  (22) 

where  Rcisth*  Fourier  cosine  transform  of  the  response, 

Rs  is  the  Fourier  sine  transform  of  the  response, 

Ic  is  the  Fourier  cosine  transform  of  the  computer  developed  response  to  impulse, 

Is  is  the  Fourier  sine  transform  of  the  computer  developed  resonse  to  impulse. 

Note  this  is  impulsive  response  at  P  due  to  unit  impulse  at  K  when  R  is  response  at  P  and  the  impact 
force  is  at  K. 


Solving  for  the  Fourier  transforms  Fc  and  Fs  yields 

4"  ^e/c 

Fc  * - i — n —  —  /r(w),  Fourier  cosine  transform 

c  /<?  +  /i 


FSU u);  Fourier  sine  transform. 


(23) 

(24) 


For  a  forcing  function  fit)  undefined  for  t  <  0  one  need  not  consider  fit)  before  t  «*  0  and  may 
choose  fit)  to  be  even,  odd,  or  neither.  The  impact  force  may  then  be  obtained  from  Eqs.  (11),  (15) 
or  (18).  Typically  one  would  compute  both  Fs(w)  and  Fc(to)  and  depending  on  which  transform  con* 
verged  toward  zero  faster  choose  either  the  inverse  cosine  transform  (Eq.  05))  or  inverse  sine 
..ansform  (Eq.  (18))  to  get  the  time  history  of  the  impact  force. 

The  process  is  then  repeated  for  the  zth  point,  etc.,  and  a  set  of  estimates  can  be  built  up.  The 
spread  present  in  them  is  then  a  measure  of  how  well  the  computer  model  of  the  structure-equipment 
fits  the  real  world. 

COMPUTATION  OF  FOURIER  SINE  AND  COSINE  TRANSFORMS  AND  INVERSE 
TRANSFORMS 

In  Ref.  (1]  it  is  shown  that  the  Fourier  cosine  and  sine  transforms  computed  over  a  finite  time  f0 
of  a  function  f(t )  expressed  as: 

Fc(m)  m  fj*  fit)  cos  mt  dt  (15) 

Fs(m)~  fj*  At)  sinvtdt  (W 

may  be  computed  from: 

Fc(m)  -  -  *(  7»)  cos  w  T,  ♦  -  ■  t&Li!.  #n*r*  ^  >  ^ 

Fs(m)  m  -  cos  mT9- X(T*)  sin  mt< 

*  •  -s--.  vr's:,--.  *  '  *  :» 
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where  X  is  the  displacement  response  of  a  linear  oscillator 
X  is  the  velocity  response  df  a  linear  oscillator 
m  is  the  circular  frequency  (rads/s)  ^ 

Tq  is  the  end  of  the  time  interval 

The  problem  of  finding  the  transforms  resolves  itself  into  finding  X(T0)  and  X(T{>)  for  a  given 
value  of  at.  T0  and  /( To)  are  known.  The  process  is  repeated  for  a  number  of  evenly  spaced  values  of 
at  over  the  desired  range.  The  method  to  find  X(T0)  and  X(T0)  is  the  numerical  integration  scheme 
presented  in  [1]  but  with  some  trigonometric  substitutions  to  prevent  computer  precision  problems. 
The  integration  scheme  is  as  follows: 

.  .  X„  sin  w'A  2S„  sin*  («k/2) 

Xg+i  -  X„  cos  m  k  +  -  -  — a. 


-  $S. 


co s*(mA/2) 

-1 


at2h 


sin  a»h 


a »JA2 


X,+i  - -atXg  sin  oth  +  X,  cos  «A  -  S„ 

-ss^ 


2  sin2  (ath/2)  sin  mh  I 

m2h2  "  2 mh  | 

where  h  is  the  time  increment,  and  S„  -  /.+l  -  /„  and  SS„_i  -  /„+,  -  2f„  +  fm.\. 

The  inverse  Fourier  cosine  and  sine  transforms,  defined  by: 

2  rwa 

frit)  m  —  I  Fr(«)  cos  dw 

n  Jo 
2  tw s 

"  — •  J0  /*(»)  sin  «r  d« 


(29) 


(30) 


(31) 

(32) 


may  be  compute^  in  an  exactly  analogous  way  by  letting  »  become  t,  T0  become  Wo,  and  h  becoming 
the  increment  in  at  in  Eqs.  (29)  and  (30)  and  multiplying  the  final  result  by  2/w. 

The  numerical  integration  method  has  two  types  of  error  present:  inherent  and  round  off.  For 
transforms  involving  functions  which  have  a  finite  number  of  finite  discontinuities  and  can  he  exactly 
described  by  a  set  of  straight  lines  or  parabolic  arcs  the  method  has  no  inherent  error.  The  numerical 
integration  equations  for  the  case  above  will  give  exactly  the  transform  value  except  for  mind  off  error 
at  any  at  regardless  of  the  increment  size  A(Af).  For  other  functions  the  closer  the  approximating 
curve  ties  to  the  true  function  the  more  exact  the  transform.  '  v‘  '  ~ 

An  alternative  method  to  compute  Fourier  transforms  is  the  fast  Fourier  transform  or  FFT.  Some 
drawbacks  are  that  with  the  FFT  the  function  must  be  long  enough  to  get  the  desired  frequency  incre¬ 
ment  Aw  and  the  number  of  frequencies  or  simples  must  be  s  power  of  2.  With  this  method  one  sim¬ 
ply  chooses  a  A<*  end  the  number  of  frequencies  desired;  or  one  may  compute  the  transform  at  any 
desired  frequency.  With  the  FFT  the  frequencies  are  restricted  to  multiples  of  A«  end  one  cannot 
obtain  transform  values  at  frequencies  in  between  samples  except  by  decreasing  the  sample  spacing  or 
interpolation. 


TWO  DEG1ES  or  rftXZDOtf  SYSTEM  TO  TEST  FOftCES  StONATUEE  ttCO*8TtUCT!ON 

In  order  to  test  the  method  to  reconstruct  the  fores  signature,  a  2  degree  of  freedom  model  ss 
shown  hi  Fig.  1  consisting  of  2  coupled  maafr^prftbMuRpar  system*  Is  employed  The  second'  Cassis 
driven  by  a  force  of  known  shape  and  the  response  of  both  masses  are  calculated  to  this  force.  In  addi¬ 
tion  the  responds  of  the  masses  to  a  unit  hnpeist  is  also  calculated.  Knowing  the  response  to  a  foroe 


For  the  following  test  cases  the  following  parameters  are  used:* 


Mi  -  1  Ib-sVin  Ki  -  27000  Ibs/in  C2  -  100  Ib-s/ia 

A/,  -  2  lb-sVin  Kx  -  9000  Ibs/itt  C,  -  200  Ib-s/in 

Case  I:  Parabolic  Force 

A  parabolic  input  force  of  form:  ...  ■  >  ^ 

Fit)  -  16000  it- 40  t2)  ^  0  <  t  <  .025  s 

Fit)  -  0  t  >  .025  s  . :t  'i  s"; 

as  shown  in  Fig.  2a  is  applied.  The  exact  analytical  solution  is  obtained  for  yx  and  by  the  method  of 
Laplace  transforms.  The  solution  is  straight  forward  but  lengthy  and  will  not  be  shown  here. 

The  solution  for  y\  and  y2  where  the  force  is  a  unit  impulse  is  also  obtained  by  Laplace 
transforms.  Using  a  computer  program,  the  sine  and  cosine  transforms  of  the  response  of  yt  to 
impulse  and  response  to  the  parabolic  force  are  computed  from  Eos.  (27*30)  where  the  endef  the  dme 
interval  f9  is  .3  s,  h  the  time  increment  is  .0025  s  and  the  frequency  m  ranges  from  0  to  2000  nda/s  in 

increments  of  S  rads/s.  These  transforms  are  then  combined  according  to  Eqs.  (23)  and  <24)  ao«0tai* 

the  sine  and  cosine  transforms  of  the  force  shown  in  Figs.  2b  and  2c.  The  inverse  transform*  are 
obtained  using  the  same  equations  as  for  the  forward  transforms  Eqs.  (27-30)  where  the  end  of  die  fre¬ 
quency  interval  W0  is  2000  rads/s,  Aw  is  5  rad/s  and  the  dme  ranges  from  0  to  .1  s  in  iMtiUMtoftoftiV 
.0025  s.  The  inverse  transforms  are  combined  according  to  Eq.  (11)  to  obtain  the  force  dme  history 
shown  in  Fig.  2d.  By  comparing  Figs.  2a  and  2d,  it  can  be  seen  dim  the  input  feme  and  the  recon¬ 
structed  force  agree  very  well.  The  same  agreement  hi  obtained  if  the  response  of  y2  is  used  toeeoee* 
street  the  force.  It  should  be  noted  that  in  order  to  reconjtniet  die  force  out  to  time  rit  ia  nenemarywr 
have  the  responses  out  to  a  latter  time  t  +  71  For  example,  in  the  case  above  the  impulse  carpels  and 
response  to  force  war e  used  out  to  .3  s  to  reconstruct  the  force  time  history  out  to  .1  s. 

Case  II:  Hewn  Ramp  Farce 

A  down  ramp  input  force  of  form:  s 


Fit)  -  100.  -  (100/ .025)  r  0  <  i  <  .025  s 
Fit)  -  0  t  >  .025 1 

as  shown  in  Fig.  3a  is  applied.  In  order  to  avoid  recalculating  long  analytical  solutions  for  each  input 
force  a  fourth  order  Rungt-Rutta  numerical  method  is  used  to  compute  the  response  tw  the  above 
force;  however  the  analytics!  solution  for  die  impulse  respowm  is  used.  Transforms  of  the  same  length 
and  increment  tire  as  In  Case  1  are  taken  end  dm  mme  procedure  is  used  to  obtain  the  font  tmtofoms 
shown  in  Fife  3b  and  3c.  The  cosine  transform  «  expected  convenpss  much  factor  than  tin  sine 
tnutsform.  3^  dm  fgptb  t^^  $$$&&  recongtoere  dgt. 

IwCi.  nfttffv  Jm  MOWS  M  HfTOfvQ  OMNI  By  U0HS  MSf  MS  C0WS  MMl  » IMMSCl  IM 

input  force.  TWs  bnprovsmem  • 9  ushtg  only  die  cosine  transform  is  for  dw  foiowtof  reatwt.  Ifdto 
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reflection  of  tile  down-ramp  function  around  the  y  axis  is  considered,  the  function  becotnesevea  and  so 
only  the  cosine  transform  is  required  in  the  inverse  transform.  When  the  inverse  transform  is  taken  by 
using  Eq.  (1$),  the  cosine  transform,  because  it  converges  faster  than  the  sine  transform  reconstructs 
the  function  for  t  >  0  more  accurately  than  both  transforms  together  as  in  Eq.  11.  Since  the  forcing 
function  doesn’t  exist  before  r  *  0  the  inverse  transform  should  not  be  considered  before  t »  0;  in  the 
above  case  the  ramp  before  r  «•  0.  For  an  unknown  force  both  the  sine  and  cosine  transforms  may  be 
examined  and  the  one  which  converges  faster  may  be  used  for  the  inverse  transform  for  t  >  0. 

Case  III:  Damped  Sloe  Force  — 

A  damped  sine  wave  input  force  of  form: 

Fit)  -  (1.  -  t/T)  sin  4 wt/T  0  <  t  <  .025  s,  T  -  .025  s 
Fit)  -  0  t  >  .025  s 

as  shown  in  Fig.  4a  is  applied.  As  in  Case  II  the  response  to  the  force  is  calculated  using  a  Runge- 
Kutta  method.  The  same  procedure  as  with  Case  II  is  used  to  calculate  the  sine  and  cosine  transform 
of  the  force,  shown  in  Figs.  4b  and  4c.  The  force  reconstructed  by  using  both  transforms  is  shown  in 
Fig.  4d. 

Case  IV:  Exponential  Decaying  Force 

A  decaying  exponential  input  force  of  form: 

Fit )  -  100  e'500'  0  <  t  <  oo 

as  shown  in  Fig.  5a  is  applied.  As  in  Case  III  the  response  to  the  force  is  calculated  using  a  Runge- 
Kutta  method.  The  same  procedure  as  in  Case  III  is  used  to  calculate  the  sine  and  cosine  transform  of 
the  force  shown  in  Figs.  5b  and  5c.  The  force  reconstructed  by  using  both  transforms  is  shown  in  Fig. 
5d.  The  force  reconstructed  by  the  cosine  transform  above  shown  in  Fig.  5e  is  closer  to  the  input  force 
because  the  cosine  transform  converges  faster  than  the  sine  transform. 

SUMMARY 

An  accurate  method  has  been  presented  to  compute  the  force  signature  of  a  body  impacting  on  a 
linear  elastic  structure.  The  various  input  forces  applied  to  a  two  degree  of  freedom  system  were  all 
reconstructed  with  good  accuracy,  and  assuming  that  that  accurate  impulse  responses  could  be  obtained 
the  method  would  work  equally  as  well  on  a  complex  multi-degree  of  freedom  structure. 
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